Spinon deconfinement above a finite energy gap in 
2d Quantum Heisenberg Antiferromagnets 
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The familiar spin- 77 quantum Heisenberg antiferromanget in a 2d square lattice is shown, within 
the non linear sigma model approximations, to be another novel state of matter that has excitations 
with fractional quantum numbers above a finite energy gap. The 1-skyrmion with an energy m 2ttJ 
is shown to be made of two 'deconfined spinons' or 'SU(2) vortices'. The many skyrmion operator 
and the wave functions that we have found are strikingly similar to quantum Hall quasi particle 
operators and wave functions. We also predict the presence of finite energy 'spin-S spinon' for a 
general spin-S Heisenberg antiferromagnets in 2d. Some consequences are briefly discussed. 



Neutral spin-i fermionic excitations now called 
'spinons' Q were conjectured to be present in the quan- 
tum spin liquid states of spin-half 2d Heisenberg antifer- 
romagnets, by Anderson in 1987. While this was readily 
shown to be present in the RVB mean field analysis of 
the spin liquid vacua and short range RVB states j| , their 
presence as a finite energy deconfined excitation in the 
ordered antifcrromagnetic vacua has remained unclear. 

A conjecture by Dzhyaloshinskii, Wiegman and 
Polyakov || and the early work of Wilczek and Zee || 
suggested statistics transmutation of skyrmion in ordered 
spin-i Heisenberg antiferromagnets in 2d. Anderson, 
John, Doucot, Liang and the present author Q however 
conjectured that finite energy 'half-skyrmion' or 'meron' 
are the 'deconfined' spinons based on some heuristic 
arguments. 

Recent inelastic neutron scattering results covering 
a large energy and momentum range in the insulating 
cuprates ||, and also Raman |jj and infrared measure- 
ments also makes the search for any signatures of spinons 
at low and high energies very meaningful and urgent. 

The aim of the present letter is to study the spectrum 
of quantum Heisenberg antiferromagnet in 2d and look 
for deconfined spinons above a finite energy gap, within 
the 0(3) non linear sigma (NLS) model approach. We 
look carefully at Belavin Polyakov n-skyrmion static 
classical solutions. The mathematical structure of the 
n-skyrmion solution in one particular parameterization 
readily suggests that each skyrmion is made of 2n 'con- 
stituent point particles [Qjl^Jl^]. 

To understand the quantum dynamics of these con- 
stituent particles we construct the skyrmion operator for 
our spin-i Heisenberg antiferromagnet and discover that 
the creation operator for these constituent particles have 
mathematical structure similar to Laughlin's quasi hole 
and quasi particle operators of quantum Hall effect. By 
a Berry phase analysis we show their spin to be half. 
Asymptotic form of the modulus of the n-skyrmion wave 
function in terms of collective coordinates is also found. 

Our starting point is the spin-S quantum Heisenberg 
antiferromagnet in a 2d square lattice with nearest neigh- 



bor interactions. Following the standard derivation IrJ 
one arrives at the 0(3) NLS model action along with the 
important lattice sum of Berry phases: 



S = 



f dxdy dt([d^n{v)] 2 - -l[c) t n(r)] 2 ) +iS B [n] 



(1) 



Here fi = x,y and n(r) is a normalized (n(r) ■ n(r) = 1) 
sub lattice magnetization vector. 

The coefficient p Q w J for S = \ case, and v s is the 
spin wave velocity. The Berry phase term is a lattice sum 



S B [n] = 27rS^(-l)l m l+l"lfl (rw) 



(2) 



Here the integers (m,n) stand for a lattice site and 
^°( n m,«) = / *^n m .n ■ {d t n m n x d u n m n ) is the single 
site Berry phase. 

Inspired by the conjecture of reference (6) we look at 
the finite energy topological solutions in our search for 
spinons. The multi skyrmion solutions found by Belavin 
and Polyakov |l(j] are extended objects with non trivial 
topology. Their Eucledian solutions of the 0(3) model in 
(1 + 1) d become the time independent classical solutions 
of our 0(3) model in (2 + 1) d. 

An n-skyrmion solution is given by 



n 

»M=ii(- 



(3) 



Here z = x + iy. The n complex co-ordinates a i 
and b i characterize the skyrmion solution. The func- 
tion w(z) and the sub lattice magnetization n(r) = 
(sin <p(r) cos 0(r), sin </>(r) sin 0(r), cos </>(r)) are related by 



w(z) = cot ■ 



(r) 



(4) 



The n-antiskyrmion solutions are obtained by replac- 
ing z by z. In our convention the spins at infin- 
ity n(oo) = (1,0,0) for the ground as well as the 
skyrmion/ antiskyrmion states, since w(oo) = 1. 
The energy of the n-skyrmion solution is given by 



1 



(5) 



a constant independent of the skyrmion coordinates 
{a i ^b i \ ! This means that at the classical level the n- 
skyrmions do not interact; thus the set {a il bj} represents 
the 2n 'zero mode' (2 dimensional) co-ordinates of the n- 
skyrmion. We will call the a and b coordinates of the 2n 
constituent particles of an n-skyrmion. 

Physically the co-ordinates a and b represent the 'cen- 
ters' of local vortex distortions of the xy compounts of 
the 3 vector field n. However, unlike the ± vortex pair 
of an xy model, there is a freedom in the choice of the 
positions of the vortex pair; for a 1-skyrmion any two di- 
ametrically opposite points on the circle \z— s 4^- \ = \ \ 
can be chosen as centers of the vortex pair, with the lo- 
cal vortex distortion occurring in appropriately rotated 
2d plane in the spin space. It is remarkable that the 
extended spin twists of a 1-skyrmion, which is more of 
a'ring like' distortion of the n field, than two point vor- 
tices, is characterized by just two points a and b, which 
we will later elevate to the level of particle degrees of 
freedom in the plane. We say that a n-skyrmion is made 
of 2n 'SU(2) vortices'. 

It should be remarked that the Belavin Polyakov's ex- 
act solutions describe either skyrmions or antiskyrmions. 
No exact solution containing both skyrmions and anti- 
skyrmions exist. In approximate solutions containing 1- 
skyrmion and 1-antiskyrmion there is a dipole like inter- 
action between them even at the classical level 0. 

The two SU(2) vortices of a 1-skyrmion state are finite 
energy solutions and they have twisted their way out of 
the log \a — b\ energy dependence, characteristic of xy vor- 
tices in 2d. We will see later that quantum fluctuations 
produce repulsion between the two SU(2) vortices of a 1- 
skyrmion ! Further, once an xy anisotropy is introduced 
the SU(2) vortices degenerate to two xy vortices with an 
energy diverging as Jm j Jz log \a — b\ as \a — b\ — > oo. 

An n-skyrmion/antiskyrmion carries a topological 
'quantum number' or a 'winding number' ± n, the degree 
of the map (x,y) — > S 2 of the field n(x,y). Here S 2 is 
the order parameter space. The 'topological density' q(r) 
is given by 



<?( r ) = ^ n ( r ) ' (^ n ( r ) x d y n ( r )) 

{d z w(z)dzw(z) — d z w(z)dzw{z)) 
= ir{l + \w{z)\ 2 ) 2 



(6) 



And J dxdy q(r) = n is the winding number. 

To keep comparison with other formalisms, the above 
topological density is related to the U(l) magnetic flux 
of the RVB gauge field fTJ: 

Si ■ (Sj x S fc ) = i(TijT jk Tki - T ik T k] Tji) ~ e l f A ™ b ' dl (7) 
where T i} = J2 a 4* c ja and n n + "a = 1 



To understand the quantum dynamics of the skyrmion 
excitation and its constituent particles for our spin-i 
Heisenberg model on the square lattice, we first construct 
the operator for the 1-skyrmion on the lattice. (No ex- 
act classical solution is available for the 'lattice version' 
of the NLS model in 2d; when \a — b\ >> 1, in units of 
the lattice parameter, the continuum skyrmion solution 
should be a good approximation to the lattice solution). 
The spin wave ground state of our spin-^ Heisenberg an- 
tiferromagnet is 



j]~(-l)"*exp{ 



i<j 



(8) 



where zi = rrii + i Ui stand for the lattice co-ordinates of 
the down spin electrons in the square lattice contain- 
ing N sites. This is a hard core boson (or equivalently 
s z ) representation of the spin wave ground state, where 
the Jastrow factor f(\zi — Zj\) ~ \ z -z-\ ar i se from the 
zero point quantum fluctuations of the gapless spin wave 
modes of the ordered antiferromagnet. And (— 1) Na is 
the Marshall sign factor, where N A counts the number of 
down spins on sub lattice A. The ODLRO exhibited by 
the above ^ hard core boson wave function represents 
the long range antiferromagnetic order, with a finite sub 
lattice magnetization m, lying in the xy plane. We can 
fix the direction of m in the xy plane by giving an addi- 
tional coherent superposition with respect to S^ otal (the 
total number of hard core bosons) in the above wave 
function Jl5| ]. 

After some algebra the co-ordinate representation of 
the 1-skyrmion operator is shown to be 



68m>-ii(!^) 

i—1 



(9) 



While w(z) of equation (3) represents the distortion of 
the sub lattice magnetization in space, equation (9) pro- 
vides the operator that produces the distortion. To prove 
the above, it is convenient to put the sub lattice magneti- 
zation along the + x-axis in the ground state and use the 



>r representation, Og}(a, b) 



^(a ) &)=nSi(UEf)'" (Note, 
unlike equation (9), here z^s run over all sites !). 

This is an important result, which shows that a 1- 
skyrmion operator actually creates one Laughlin 
1 16] quasi- hole and one quasi-electron like object, 
when acting on the hard core bose fluid. Similarly an n- 
skyrmion state can be created by a product of n of Laugh- 
lin quasi-holes and n of quasi electron like operators. The 
operator for anti skyrmion is obtained by replacing {z^ 
by {z^ in equation (9). Our quasi hole operator exactly 
coincides with Kalmayer-Laughlin's spinon operator [pi. 

Since at the classical level, the a and b co-ordinates 
of the n-skyrmion solution represent the zero-mode co- 
ordinates, it is legitimate to elevate them to the level of 



2 



new low energy collective co-ordinates or 'particle degrees 
of freedom as in the case of Laughlin's quasi holes 
and quasi electrons. 

While the a and b 'quasi particles' (SU(2) vortices) do 
not have any interaction at a classical level, the modified 
quantum fluctuations in the the presence of a skyrmion 
can induce interactions among the constituent particles. 
The induced interaction V qu (a, b) is obtained from the 
difference in zero point energy of the spin wave modes in 
the presence Ku)^' b ^ and absence hco^ of a 1-skyrmion: 

V qu (a, &)^£ -\H h ^ (10) 

Using the method of Rodriguez Jl?]] and also Marino [jl8| , 
who did not use the constituent particle interpretation, 
we have calculated the induced interaction between an a 
and b particle. We find that the a and b particles always 
repel. It is energetically advantageous for the constituent 
particles to be infinitely apart and reduce the zero point 
fluctuation energy. In particular the quantum fluctuation 
corrected energy of a 1-skyrmion state is: 

^= 1 (|a-6Hoo)«2 7 rJ 

This proves that the two constituent particles are indeed 
deconfincd. 

I have obtained an asymptotically exact modulus of 
the wave function of the 'ground state' of the n-skyrmion 
state, for large separation of the a and b coordinates as 

i*„ [o, b] i ~ n \ ai - a > 1 1 6 * - b > i n - b i i ^ ( n ) 

The method I have devised involves expressing the path 
integral for vacuum to vacuum amplitude in a way that 
gives the desired wave function in terms of the Jaco- 
bian of transformation from the function w(z) to the 
zero mode coordinates. A plasma analogy shows that 
the a and b particles are indeed unbound, giving another 
proof for the dcconfmcmcnt of the a and b particles in 
the skyrmion state. 

The spin of a skyrmion or the constituent SU(2) vor- 
tices and their exchange statistics are rather ill defined 
in view of the large quantum fluctuation contained in 
the spin wave ground state. The 'mathematically sharp' 
point particles a and b are extended physical objects with 
power law form factors. (This is similar to the case |l9| of 
vortices in 2d superfluid He 4) . The integrated missing or 
excess density around a vortex gives the z-component of 
the spin of the a and b particles. This calculation turns 
out to be hard with our hard core bose fluid. So we take 
a Berry phase approach. 

In our Berry phase analysis the spin of the con- 
stituent particle or the SU(2) vortices appears as a 'chain 



anomaly'. We consider a one skyrmion solution and an- 
alyze the lattice sum of the Berry phases under a global 
rotation of all the spins through 2tt about x-axis, the 
direction of sub lattice magnetization at infinity. If the 
spin projections of the constituent particle along the x- 
axis are a a and <Jb, the corresponding lattice sum of Berry 
phases should contribute 2iro~ a and 2nab for the a and b 
particles in a spatially separated fashion. 




field n(x,n;t) for an/0 chain. The 'hole' just vanishes 
for the n = chain; inset shows its n(x, 0; t) at t = 0. 

We find an interesting phenomenon which we call a 
'chain anomaly'. We group, following Haldane pT3| ] , the 
square lattice sum of Berry phases into sums over chains, 
running parallel to the x-axis: 

S B = 2nSj2(-^l d (12) 

n 

where the n-th chain Berry phase is given by the sum 

nlf = ]T (-) m o°(n m ,„) 

m=-JV,.JV-l 

= J -^- n (x,n;t) ■ (d x ii(x,n;t) x d t n{x , n; t)) (13) 

The chain Berry phase counts the winding number of 
map (x, t) —¥ S 2 of the field n(x, n; t) for n-th chain. 

We choose the length of our lattice along the x-axis 
to be 2N, where N is an odd integer. For convenience 
we choose (0, 0) to be the center of the square lattice 
and the two skyrmion co-ordinates to be the lattice sites, 
a = (to, o) and b = (—to, 0) on the x-axis at time t = 
and T. Our global rotation amounts to giving the follow- 
ing time dependence a(t) = me'~ and b(t) = —me l ^r~ 
to the constituent particle co-ordinates. 

During the time evolution, the spin field n(x, n; t) of 
the 1-skyrmion state of any chain n attempts to wrap the 
S 2 sphere. All, except the x-axis chain ( n = 0) fail to 
wrap the S 2 sphere completely. They all leave a hole, 
as shown in figure 1 , making the winding number (Berry 
phase) identically zero. The x-axis chain just manages to 
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wrap the S 2 sphere, once for the chain sum —N, .. — 1 and 
again in the opposite sense for the chain sum 0, 1, ..N — 1, 
thereby contributing a phase of it for the left half and 
—tt for the right half of the chain. We identify these two 
phases, which arise predominantly from the a-vortex re- 
gion and b- vortex region, with the spin-i Berry phases 
of a and b particles. 

The lattice sum of Berry phases can be grouped in 
many different ways and we always find one singular 
chain passing through two diametrically opposite points 
of the circle \z — = of our 1-skyrmion solu- 

tion. The two ends of the chain either close on themselves 
or go off to infinity. The spin field n for these singular 
chains trace a closed curve in S 2 connecting two antipo- 
dal points at any given time; it undergoes 'ir-twist' twice, 
corresponding to the presence of two spinons. For exam- 
ple, the circle \z — 9 4j^-\ = \ 3 ^\ itself is a singular chain, 
for which the n field traces a great circle on S 2 at any 
given time. 

In spin-i Heisenberg chain each 7r — twist of n field in 
space corresponds to one spinon p3|. In a sense, in the 
1-skyrmion configuration a 'Id-chain' or a 'string' con- 
taining two spinons is embedded in the 2d plane in an 
irreducible and dynamical fashion. 

Our chain anomaly is missed if we use Haldane's argu- 
ment |H| for the calculation of the Berry phase. Haldane 
suggested that if the n(x, y; t) field is continuous and non 
singular (as it is for the case for n-skyrmion solution), the 
chain Berry phases should be all identical and as a result 
the staggered sum should be identically zero, for an even 
number of chains. However, we find that our singular 
chain containing the a and b co-ordinate is an exception 
to this and it contributes a Berry phase of it and — tt to 
the two particles a and b. 

The two SU(2) vortices or the two spinons of a 1- 
skyrmion state, by the argument given above, carry spin 
half projections of value 5 and — \ along the x-axis. In 
addition the spinons carry +1 chirality quantum number 
because they are constituent particles of a + 1-skyrmion 
state. Thus a spinon in addition to a spin quantum num- 
ber carries a ±1 chirality quantum number. As the anti- 
ferromagnetically ordered ground state has zero chirality, 
the finite energy skyrmion states are chiral doublets. It 
should be interesting to look for this degeneracy in the 
excited states in numerical studies of finite systems. This 
degeneracy is also reminiscent of the degeneracy of the 
spinon states in the 1-d Haldane-Shastry model arising 
from the Yangian symmetry. 

Our analysis goes through for any spin. In general for 
the spin-S square lattice Heisenberg model we get decon- 
fined 'spin-S' spinons or the SU(2) vortices. 

The spinon deconfincment that we have found has sev- 
eral interesting consequences which we hope to discuss 



in the future: i) suggests a way to produce charge 2e 
skyrmions to get superconductivity in the doped Mott 
insulator, ii) as the energy of a deconfined spinon w n J, 
its possible signature in the top of the spin wave band as 
well as in the infrared and two magnon Raman measure- 
ments, iii) short range RVB state viewed as the conden- 
sation of skyrmion and antiskyrmion in the ground state; 
chiral symmetry broken Kalmayer-Laughlin like states as 
condensation of unequal density of skyrmion and anti- 
skyrmion and iv) consequence of our present picture to 
skyrmion doping in quantum Hall ferromagnets. 

I wish to thank R. Shankar for discussions and M.B. 
Silva Neto for bringing reference (17) to my attention. 
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